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The thick center vortex model with the idea of using domain structures is used to calculate the 
- - - potentials between two G(2) heavy sources in the fundamental, the adjoint and the 27 dimensional 

(N- representations. The potentials are screened at large distances. This behavior is expected from 

, the thick center vortex model since G(2) has only a trivial center element which makes no 

>— ^ ■ contribution to the average Wilson loop at the large distance regime. A linear potential is obtained 

at intermediate distances for all representations. This behavior can be explained by the thickness 
{3JQ[ of the vortices (domains) and by defining a flux for the trivial center element of G(2). The role of 

^ ■ the SU(3) subgroup of G(2) in the linear regime is also discussed. The string tensions are in rough 

agreement with the Casimir operators of the corresponding representations. 
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I. INTRODUCTION 



The last candidate for the fundamental theory of hadronic force is the quark model. The interaction between quarks 
(j^ ■ is described by the exchange of the non-Abelian gauge fields called gluons. Quarks were introduced by Gellman and 
. Zweig for the first time in 1964. Besides the successfulness of the quark model, no free quark has been observed 
yet. However, the particles in the Lagrangian of any theory must exist in the physical spectrum. Now that the 
. quarks appear in the QCD Lagrangian, the main question is as follows: Where are the quarks? Are they permanently 
^ ■ confined? In fact, these are not the quarks which are not found in the nature but it is the color which is confined 
. and quarks are confined because of their colors. In other words, particles with color degrees of freedom are confined. 
\G ' Therefore, all free particles are colorless, and colored particles cannot be found in the particle spectrum. Also gluons, 
particles in the adjoint representation of the gauge group of QCD have colors, and are absent in the particle spectrum, 

■ as well. 

04 ' Let us define confinement from the point of view of the potential in a quark-antiquark pair. At short distances the 

■ potential in a quark-antiquark pair is Coulombic and can be explained by the gluon exchange between quarks. Since 
] ' the coupling constant is small, perturbative methods work very well in this region. As the distances between quarks 

" . [ , increase, the coupling constant increases as well and perturbative techniques do not work anymore. This happens 
^ ■ at an intermediate regime where the potential between quarks rises linearly until the energy reaches a point where 
\^ [ it is large enough to create a quark-antiquark pair in the vacuum. Then, these secondary pairs couple to the initial 

■ ones and screen them. However, by removing the dynamical quarks from the theory, the creation of the secondary 
" quarks does not happen and the potential between the initial quarks rises linearly with distance. Figure ([1]) shows 

a typical behavior of the potential between static quark s, a Coulombic potential at short distances, and a linear 
potential for intermediate and large distances. For this figure, dynamical quarks are removed from the theory. The 
linear behavior is confirmed by Regge trajectories in the particle spectrum [l[. The flux tube created at this linear 
regime has nonlinear nature, and it is not possible to use ordinary perturbation methods to study its characteristics. 
Therefore, to investigate the quark confinement phenomena, one should use nonperturbative methods. Lattice gauge 
theory has been very successful to reproduce the potential between pairs of colored particles like quarks and gluons. 
Other approaches used to study confinement include phenomenological models. In these models, the QCD vacuum 
is filled with some topological configurations confining colored objects. The most popular candidates among these 
topological fields are monopoles and vortices. Other candidates include merons, calorons,eic. There are very strong 
correlations between these various objects, though. Therefore, the mechanism of confinement has not been understood 
by a unique procedure yet, and it is one of the unsolved mysteries of quantum chromo dynamics. 

In this paper, we study the center vortex model within the G(2) gauge group. The center vortex model was initially 
introduced by 't Hooft [2|. It was able to explain the confinement of quark pairs, but it was not able to explain the 
intermediate linear potential, especially for higher representations. The model has been improved to the thick center 
vortex model by M. Faber et al. Even though it is a simple model, the thick center vortex model reproduces 
the intermediate linear potential qualitatively, in agreement with Casimir scaling and the asymptotic large distance 
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potential which is proportional to the A^-ality of the gauge group. For intermediate distances, one expects to see a 
linear potential proportional to the distance between quarks, in agreement with lattice calculations. Based on lattice 
results [3, 0] , the coefficient of the linear part which is called the string tension, is proportional to the eigenvalue of 
t he quadratic Casimir operator of the corresponding representation. This is called Casimir scaling. The Casimir 
scaling regime is expected to extend roughly from the onset of confinement to the onset of screening. The A'^-ality 
of each representation is given by mod(ri — m), where n is the number of quarks and m is the number of antiquarks 
constructing the representation. At large distances, where the distance between the quark and antiquark increases, 
the energy stored in the string increases and gluon pairs are created from the vacuum. These gluons screen the initial 
sources to the lowest dimension with the same iV-ality. In this regime, sources with the same A^-ality obtain the same 
string tensions. Those with zero A^-ality are completely screened. 

In the vortex model the QCD vacuum is filled with linelike (in three dimensions) or surfacelike (in four dimensions) 
objects, which carry magnetic fiux quantized in terms of the center elements of the gauge group. At large distances 
between color sources, the interactions between Wilson loops and center vortices lead to a linear potential proportional 
to the iV-ality of the given representation. 

Besides the successes of the center vortex model, the relation between confinement and center symmetry has not 
been understood very well,yet. To study the role of center elements, one of the attractive methods is studying 
confinement in gauge groups without nontrivial center elements, like G(2). Lattice calculations in G(2) [6] show 
color screening for large distances and a linear potential for intermediate distances. Since the center is trivial, the 
center vortex model predicts color screening. It is not so clear why a linear potential can be expected at intermediate 
distances. Because of these two features confinement at intermediate distances and screening at large distances G(2) 
gauge theory resembles SU(A^) gauge theories with dynamical quarks. These screen static quarks at large distances. 
Therefore, G(2) is considered as a mathematical laboratory to examine the properties of SU(A^) gauge theories, even 
the ones which can not be understood in SU(A^) gauge theories themselves. 

In the next section, we give a brief review of the thin and thick center vortex models. In Sec. Ill, some general 
features of the G(2) gauge group are discussed, and in Sec. IV we apply the improved thick center vortex model, called 
the domain structure model, to the G(2) gauge group. We study the possible reasons why confinement is observed in 
the G(2) gauge group based on its SU(3) subgroups. Potentials between static sources of higher representations, 14 
and 27, are calculated in Sec.V, and Casimir scaling is discussed. We end the paper with a conclusion. 



V(r) 




FIG. 1: Potential between static quarks. At small distances the coupling constant is small enough to use perturbative methods. 
A Coulombic potential is obtained in this regime. At large distances the coupling constant is large, and nonperturbative effects 
lead to a linear potential. 
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II. THICK CENTER VORTEX MODEL 



The idea of describing confinenient using vortices goes back to 't Hooft 2] . In four dimensions vortices are topofogical 
field configurations of the vacuum forming closed surfaces which can be linked to Wilson loops. Wilson loops are 
gauge-invariant observables obtained from the holonomy of the gauge connection around given loops. Confinement is 
obtained from random fluctuations in the linking number. A vortex piercing a Wilson loop contributes with a center 
element Z somewhere between the group elements U of the gauge group 

W{C) = TipUU ...U] — > TrpUU ...{Z)...U]. (2.1) 



Because center elements commute with all the members of the group, the location of Z in Eq. (|2.ip is arbitrary. The 
effect of piercings on Wilson loops can be well understood in the gauge group SU(2), which has — / as the only 
nontrivial center element. / is the 2x2 unit matrix. Vortices piercing the loop an even number of times or not 
piercing it at all, do not have any effect on W{C). Odd linking numbers change the sign of W{C). One can derive 
a formula for the string tension a assuming that vortices are thin and pierce Wilson loops in single plaquettes with 
independent probabilities /. We get 

{W{C)) = Y[{{1 -f) + .f{-l)}{Wo{C)) = cxp[-a(c)A](M^o(C)) (2.2) 

assuming that (Wo(C)) is the expectation value of the loop when no vortex pierces this loop and the string tension is 

a=-iln(l-2/). (2.3) 

The vortex model works very well for the fundamental representation, and it gives the asymptotic string tension 
correctly. It also works well for the adjoint representation at large distances, where one expects a screened potential. 
Since we have a trivial center element for the adjoint representation, the string tension must be zero according to 
Eq. ()2.ip . In other words, Z is unity and it does not change the Wilson loop. However, lattice calculations show 
an intermediate string tension for higher representations 0, S which can not be predicted by the thin center 
vortex model or Ea. (l2.ip . M. Faber et al. have improved the model to the thick center vortex model to be able to 
reproduce the intermediate string tensions. They have given a thickness to the 't Hooft vortex. Mathematically, they 
have replaced the center element Z by an element of the gauge group, G, 

W{C) TvpUU ...U] — > TrpUU ...G...U\ (2.4) 

where G is 

G{x,s) = eyip{iac{x)li ■ L). (2-5) 

Li's are the generators of the group in the representation j, it is a, unit vector, and S is an SU(A'^) element of the 
group in the representation j. adx) gives the profile of the vortex, and it depends on that fraction of the vortex which 
is pierced by the loop. It depends on the shape of the loop C and the position of the center of the vortex relative to 
the perimeter of the loop. 

If the Wilson loop is linked by m vortices, centered at positions xi, X2, a^m, then 

W{C) W[C; {x„ Si\] = TrpUU . . . UG{xa, Sa}U . . . UG{xb, Sb}U . . . UG{x^, S^}U]. (2.6) 

To obtain the average Wilson loop, W{C) must be averaged over different gauge group orientations Si. For simplicity, 
we first look at the SU(2) case. Equation (|2.5p for SU(2) is 

G{x, s) — ex-p{iac{x)~rt ■ L) = S exp{iac{x)L3)S'^ . (2-7) 

L3 is the element of the Cartan subalgebra. Independently averaging every G{xa, Sa} over orientations in the group 
manifold specified by Sa, we get 

G(a)- Jsexp[iaL3]SUS = gj[a,]l2j+i (2.8) 



where 

(7„- \fy„] = 

2.7 



9Aac] = j^TT{exp[iaLs\). (2.9) 
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Thus, averaging the Wilson loop over all Sa leads to [s*] 

(WiC)) = -f) + fg,[a,]}{WoiC)) (2.10) 

X 

= exp[^ln{(l - /) + fg,[ac]}]{Wo{C)) (2.11) 

X 

= exp[-a,A](Wo(C)) (2.12) 

where CTc is 

'^c = ^ EM(1 - ./) + fgd^c]}. (2.13) 

X 

For large loops, the center element is located completely inside the loop; thus 5j[Qc] = —1, and the string tension is 

a, = ln(l -/-/) = - ln(l ~ 2/). (2.14) 

This string tension is the same for all half integer representations j. For integer representations for which a flat 
potential is expected at large distances, the center element is trivial and gj[Q!c] = 1; thus cTc = according to 

At intermediate distances where the loop is not large enough to contain the whole vortex, Cc is 0] 

^c = {a^j(j + l) (2.15) 

where 

= ^J2 "'(^)- (2.16) 



xeA' 



Thus, for small / and small loops, the string tension is proportional to the eigenvalue of the second order Casimir 
operator. 

Equation ()2.10p can be generalized to the SU(iV) gauge theory 



N 



(WiC)) = n (1 - E /"(I - 9r[d:ixm{WoiC)) (2.17) 

X 71—1 

with 

g^[d] = -^Tr(exp[m • H]). (2.18) 

X is the location of the center of the vortex, the dimension of the representation r, /„ the probability that any given 
unit is pierced by a vortex of type n, and {Hi, i = 1,2, — 1} are the generators spanning the Cartan subalgebra. 
a"(a;) shows the flux profile for the vortex of type n. We have {N — 1) types of vortices for any SU(iV) gauge theory. 
The potential and the string tension are given by Q 

ViR) = Elnjl - E fnil-Regr[aUx)])\ (2.19) 

and 

1 

ac = -jY.^n[l-Y, fnil-RegridUx)])]. (2.20) 

X n—1 

To reproduce the potential or the string tension, one has to find the Cartan subalgebra and describe the profile 
function a^ix). The flux proflle should have the following properties: It should go to zero for a Wilson loop of length 
i? as a; — > ±oo; if the center of a vortex is completely inside the loop, it will lead to a maximal multiplicative factor 
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exp(2^Hi) Q Zn {n = 1,2, ...,N — 1) corresponding to the center elements of the group; and finally as i? the 
contribution of the vortex to the loop also goes to zero , so aii{x) — >■ 0. 

Various vortex profiles dQ{x) may be chosen. M. Faber et al. used, for the SU(2) case, 

dl{x) = 7V"[1 - tanh(ay(a;) + 4)] (2-21) 

R 

where n indicates the vortex type, a, h are arbitrary constants, and y{x) is 

yx — R \R — x\^ X 

y{x) is the nearest distance of x from the timehke side of the loop. The normalization constant TV" is obtained from 
the maximum flux condition, where the loop contains the vortex completely. 



exp(id;" • H) — z„/ 



(2.23) 



with 



e Zn 



(2.24) 



and / is the unit matrix. 

The thick center vortex model has worked well for the SU(2), SU(3), and SU(4) gauge theories 0,0. A linear 
potential in rough agreement with Casimir scaling is observed for all representations, and the asymptotic string tension 
proportional to the iV-ality of the given representation has been reported. To increase the length of the linear regime 
at intermediate distances, J. Greensite et al. @ have suggested postulating a kind of domain structure in the vacuum, 
with magnetic flux in each domain quantized in units corresponding to the elements of the center subgroup, including 
the identity element. They have applied this idea to the SU(2) gauge group which has one trivial and one nontrivial 
center element. Thus, instead of one vortex, they use two kinds of domain structures corresponding to the center 
elements, the trivial and nontrivial elements. To obtain the Wilson loop, Eq. (|2.10[) should be changed to 

{WiC)) = - h - /o) + hg-jWl] + hdiiaUmoiC)). (2.25) 

X 

fi and fa are the probabilities of plaquettes that belong to domains corresponding to the nontrivial and trivial center 
elements, respectively, and define the profile function for domain structures associated with the trivial and 
nontrivial center elements. To apply the domain structure model to the SU(A^) gauge group, one has to change n in 
Eqs. ()2.17p . ()2.19p and (|2.20|) from zero to the number of nontrivial center elements, where zero corresponds to the 
trivial center element. 

G(2) lattice gauge calculations show confinement for intermediate distances Q. The center vortex model can predict 
correctly the screening of the potentials at large distances since the center of the G(2) gauge group is trivial. But, 
what accounts for the intermediate string tensions? We apply the domain structure model, which includes the trivial 
center element, to the G(2) gauge group in Sec. IV and discuss the possible reasons for the linear potential at large 
distances. But first, in Sec. Ill, we present some general features of the gauge group G(2). 



III. GENERAL PROPERTIES OF G(2) GAUGE THEORIES 

G(2) is one of the exceptional Lie groups which is its own covering group fTo|. The rank of G(2) is 2 and it has 14 
generators, two of which are simultaneously diagonal. Its fundamental representation is seven dimensional, and the 
dimension of the adjoint representation is 14. The group G(2) is real, and it is a subgroup of S0(7) with rank 3 and 
21 generators. The determinant of the 7x7 real orthogonal matrices U of the group S0(7) is 1, and 

UU'^ = 1. (3.1) 

G(2) elements satisfy a constraint called the cubic constraint. 



Tabc — TdefUdaUebUfc 



(3.2) 
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T is a totally antisymmetric tensor, and its nonzero elements are 

2^127 = 7^154 — Ties = T235 = T264 = 7374 = T^re ~ 1. (3-3) 

Because of these constraints the number of generators of G(2) is reduced to 14. 

A quark in the fundamental representation {7} of G(2) can be color screened by gluons [loj . 

{7}(g){14}®{14}®{14} = {l}®... (3.4) 

Therefore, in G(2) gauge theories, the string tension between static quarks can be broken by gluons without the 
presence of dynamical quarks, as is necessary in SU(A^) gauge theories. 

SU(3) is a subgroup of G(2). Under SU(3) subgroup transformations, the seven and 14-dimensional representations 
of G(2) decompose into the SU(3) fundamental and adjoint representations 

{7} = {3}®{3}®{1}, (3.5) 

{14} = {8} ©{3} ©{3}. (3.6) 

The second equation may be interpreted in such a way that the 14 gluons of G(2) consist of the usual eight gluons 
of SU(3) plus six additional gluons which transform like the SU(3) fundamental quark and antiquark. One of the 
differences between these six gluons and the SU(3) quarks is that the former are bosons while the latter are fermions. 
Because all G(2) representations are real, the fundamental representation is equivalent to its complex conjugate; 
therefore quarks and antiquarks are identical. 

Simultaneous diagonal generators of this representation are [lo| 

= ^(7^11 - 7^22 - P55 + 7^66), (3.7) 

77« = -^(Pii + P22 - 2F33 - P55 - Pm + 2P77) (3.8) 
\/24 



where {Pij)ai3 = SiaSjp , and a, P indicate the row and the column of the matrices, respectively. From Eqs. p.7p and 
(|3.8[) . it can be shown that these generators have SU(3) Cartan generators, and they can be written as 



7f" = — . (3.9) 




A° (a = 3,8) are the two diagonal SU(3) generators. Since SU(3) is a subgroup of the G(2) gauge group, the G(2) 
weight diagram can be obtained by a superposition of 3, 1, and 3. Therefore, in the SU(3) subgroup of G(2), the 
center elements of G(2) can be constructed from Z(3), the center of SU(3) 




Z = \ 1 . (3.10) 



73x3 is a 3 X 3 unit matrix and z e {l,e^^} are elements of Z{i). The three Z matrices of Eq. p.lOp commute 
with the eight generators of the SU(3) subgroup of G(2) but not with the remaining six generators. This is why the 
centers Z given by Eq. (|3.10l) are not the center elements of G(2). The center elements should commute with all 
the generators. Therefore, the center of G(2) should be trivial and contains only the identity. This is true for higher 
representations of G(2) as well. Since the center element of G(2) is trivial for all representations, at large distances 
the potentials between objects in all representations are flat and screening happens. Equation p.4[) shows this fact 
for the seven dimensional representation. Another consequence is that the concept of A^-ality does not apply to the 
G(2) gauge group, in contrast to the SU(iV) gauge groups. 

In the next section, we calculate the potential between two G(2) quarks using the domain structure model. 



IV. G(2) AND VACUUM DOMAIN STRUCTURES 



Numerical calculations show a linear behavior at intermediate distances and a flat potential at large distances 
for the potential between two quarks in the fundamental and higher representations of the G(2) gauge theory. Even 
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though G(2) does not have any nontrivial center element, its trivial center element can explain the flat potential 
according to the thick center vortex model. Based on this model no vortex means no string tension at large distances, 
and this is in agreement with the flat potential for the G(2) quarks at large distances. In this section we reproduce 
the potential in a G(2) quark antiquark pair using the idea of domain structures of the vacuum in the thick center 
vortex model. 

Greensite et al. have introduced domain structures to increase the length of the linear regime in the thick center 
vortex model i]. They have suggested that the domain structure idea can be applied to the G(2) gauge group which 
has only a trivial center. With this motivation, we apply the idea to calculate the G(2) potentials explicitly. In 
SU(A^), zero 7V-ality representations have zero string tensions at large distances, while they have linear potentials at 
intermediate distances. The thick center vortex model has been able to reproduce this behavior, and we have been 
motivated to use this model to reproduce G(2) potentials which show the same behavior according to the lattice results 
but with the idea of domain structures instead of center vortices. Thickening the vortices leads to the confinement for 
higher SU(iV) representations. Since G(2) has only a trivial center element, we apply the thick center vortex model 
to this group with the constraint that the total magnetic flux measured by the Wilson loop holonomy is quantized in 
terms of its trivial center element. 

In the domain structure model, the vacuum of the Yang-Mills theory is filled with field configurations called domain 
structures, corresponding to trivial and nontrivial center elements. For example, for SU(2)there exist two type of 
domains, one associated with the trivial one, /, and one associated with the nontrivial one, —I. The main difference 
between the trivial domain and other domains is the lack of a Dirac-3 volume for trivial domains. 

Equation (|2.19p can be used in the SU(A^) gauge group for the modified thick center vortex model, if n starts from 
zero instead of 1, where n = corresponds to the trivial center element. For the G(2) case with only a trivial center 
element and no other nontrivial center elements, V{R) is given by 

V{R) = ^In {1 - /o(l - RegAdUx)])] ■ (4.1) 

X 

where /o is the probabihty that any given unit is pierced by a trivial vacuum domain, gr has the same form as in Eq. 

Tim . 

gr[d] = -|-Tr(exp[id • H]). (4.2) 

We use the fiux profile 

a°(x)=iV°[l-tanh(ay(a;) + -^)]. (4.3) 

The zero index indicates the trivial domain of G(2), and the index i is the index of the Cartan generators, a and b 

are free parameters of the model. The normalization factor is obtained by using the maximum flux condition of Eq. 
(|2.23p . For G(2), the equation is changed to 

exp(ia" -H) =1. (4.4) 

Therefore, 

exp(a7'°^i?i + <°^i?2) = / = exp(27ri). (4.5) 
Using the Cartan subalgebra of G(2) from Eq. p.9p . a™"^ and a''^"'^ are obtained. 



= 0, a^'^^ = 2ttV2A. (4.6) 



Thus, ai{x) and a2{x) are 

aUx)^0, (4.7) 



2tt 
2tt = 



^/8 



/24 



/24 



«2(a;) = - tanh(ay(x) + A)]. (4.8) 

Putting everything together, the potential between two G(2) static quarks is obtained from Eqs. (|4.1I) to (|4.3p and 
is plotted in Fig. The free parameters a, b, and / are chosen to be 0.05, 4, and 0.1, respectively. The potential 
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is screened at large distances as expected from the model. At large distances the domains are contained completely 
inside the loops, and the corresponding contribution, which is the trivial center element, does not change the Wilson 
loop; therefore, the string tension is zero. At intermediate distances a linear regime is observed. For this regime, 
some fraction of the trivial domain is located inside the loop. Therefore it can make a nontrivial contribution to the 
Wilson loop and lead to a nonzero string tension. 

To study the possible reasons for the linear part, we plot Re(gr [0^2(2;)]) versus x for the Wilson loop, with various 
spatial dimensions of R. Figure ^ plots Re((7r[a2(^)]) versus a; for i? = 100. x shows the location of the center of 
the vortex. The left leg of the Wilson loop is located at zero, and the right leg is located at x = i?. The maximum of 
Re((7r) is equal to 1 and indicates either the situation where the domain does not link to the Wilson loop or the one 
where it is located completely inside the loop. For the profile we are using, the size of the domain is proportional to 
the inverse of the parameter a, and it is about 20. Therefore, when the Wilson loop spatial length is equal to 100, we 
are sure that a complete link is established between the Wilson loop and the vortex. This fact is observed from the 
figure, since Re{gr) reaches 1 when the center of the vortex is located at a; « 50. 

Tie{gr) reaches a value of 1 when the vortex is completely inside a Wilson loop. In SU(iV) the result would be 
different, and a non-trivial center element is obtained. Another interesting feature of the figure is the minimum of 
Ke{gr) which is around —0.28. In general, for the SU(iV) cases Ke{gr) varies from 1, corresponding to the trivial 
center element, to the values corresponding to the nontrivial center elements. Since G(2) does not have any nontrivial 
center element, predicting the amount of the lower limit is not clear. The minimum of lie(gr), can be interpreted 
using the SU(3) subgroup of G(2). SU(3) has three center elements (domains): e^, where n changes from zero to 2. 
Re(gr.[a2(a;)]), corresponding to these center elements, varies between 1 and —0.5 since 



Now, calculating Re(5r.[a2(x)]) of the G(2) gauge group using its SU(3) subgroup and normalizing it with the dimension 




(4.9) 
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FIG. 2: Potentials between two G(2) static quarks obtained from the domain structure model. The potential is screened at 
large distances, in agreement with the fact that the gluons of the G(2) gauge group can couple to the quarks in the fundamental 
representation and give a singlet. A linear potential is observed for the intermediate distances. 
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of the subgroup, we get the minimum of Re((7r[«2(^)]) 

mill 



1 1 / gRia)min X 

-(Tr(e-^U,) - -Tr 1 

* gk(a)min X dr 



(4.10) 



where gR{a)min is the minimum of Re(g,.) for the SU(3) subgroup, and it is equal to —0.5 as mentioned above, dr is 
the dimension of the subalgebra which is equal to 3. Thus, the minimum of Ke{gr) for the group G(2) is obtained as 



Re{gr 
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(Tr(e^"-«)^in)- -(-1.5 + 1-1.5) 



-0.28. 



(4.11) 



Some other local minimums are observed for Re(gr) in Fig. ([3]) which may be interpreted by the SU(2) subgroups 
of G(2). Reaching the minimum of —0.28 happens for the intermediate distances R as well. Figure Q shows 
Re{gr[a^{x)]) for R = 20. 

Going back to Fig. ([2]), because of the trivial center of the G(2) gauge group, at large distances a screened potential 
is expected from thin or thick center vortex models. However, the linear regime is observed as a result of using thick 
vortices, as seen for the SU(Af) higher representations with zero N-ality. On the other hand, studying the flux profile, 
'Re{gr[ct2{x)]) , the role of the SU(3) subgroup of G(2) must be considered as another possible reason to interpret the 
linear behavior. One may claim that, at intermediate distances, the SU(3) subgroup of G(2) has a dominant role, and 
this fact leads to the observation of a linear potential at this regime. 

Pepe et al. [l3l have studied the role of the SU(3) subgroup with a Higgs mechanism. They have exploited a scalar 
Higgs field in representation 7 of G(2), which breaks G(2) to SU(3). This allows for interpolation between theories 
with exceptional confinement like G(2) and theories with ordinary confinement like SU(3). An interpretation of G(2) 
confinement via its SU(3) subgroup and using the domain structure model needs more investigation, and the present 
paper may be considered as one of the starting points. 
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FIG. 3: Re(g,. [0^2(2;)]) versus x is plotted for the fundamental representation of the G(2) gauge group. It varies between 1 
and —0.28. The value of 1 indicates the situation where the Wilson loop either does not link the vortex or it links the vortex 
completely. Re((7r[a2(2;)]) = —0.28 can be explained by the SU(3) subgroup of G(2) as explained in the text. 
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In the next section, we find the potentials between heavy sources for the higher representations 14 and 27, and we 
discuss Casimir scaling of the string tension. 

V. POTENTIALS BETWEEN STATIC G(2) SOURCES OF HIGHER REPRESENTATIONS 

Sources with dimensions 14 and 27 are obtained by 



To calculate the potentials with the domain structure model, one has to use Eqs. (|4.1[) and (|4.2[) . As Eq. (|4.2p 
demands, H , the Cartan subalgebra of G(2) must be calculated for the 14- and 27- dimensional representations. We 
use the weight diagrams of the SU(3) gauge group, and we decompose each representation of the G(2) gauge group 
to its SU(3) subgroups to obtain the weight diagrams for G(2) representations. It is clear from Eq. p.9p that the 
weight diagram of the G(2) fundamental representation can be obtained by Fig. ([5]). This is because the fundamental 
representation of G(2) can be decomposed into its SU(3) subgroup representations 



The 8 indicates the adjoint representation of SU(3). Thus, the Cartan subalgebra of the adjoint representation can 
be calculated from the weight diagrams of Fig. ([6]) and Eq. ()5.3p . 



{7} ® {7} = {1} ® {7} ® {14} ® {27}. 



(5.1) 




(5.3) 



(5.2) 
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FIG. 4: The same as Fig. ((3)l but for R = 20. The minimum of Re(gir-[a2(2;)]) reaches —0.28 for this intermediate R as well. 
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FIG. 5: Weight diagram of the fundamental representation of G(2) decomposed into its SU(3) subgroups. 



where A3, A3, and Ag are the Cartan subgroups of SU(3) for the representations 3, 3, and 8, respectively. 

Representation 27 of the G(2) gauge group is obtained from Eq. ()5.1|) . But hke the fundamental and adjoint 
representations of G(2), the weight diagram of representation 27 can be obtained by decomposition to its SU(3) 
subgroups. Figure ([7]) shows this decomposition, where 6 and 6 are the six dimensional representations of SU(3). In 
fact, 



{27} = {8} ® {6} ® {6} © {3} ® {3} © {!}. 
Therefore, the Cartan subalgebra of representation 27 is 
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(5.6) 



We recall that for all representations of G(2) we use the normalization condition 

Tr[TaTb] = ^Sab- 



(5.7) 



The next step is to find the normalization factors in Eq. (|4.3p for these two new representations. We use Eq. (|4.4p . 
with H the Cartan subgroups of representations 14 and 27. For both representations a™°^ = and a™°^ is equal 
to 7rA/24 for the adjoint and 67r\/24 for the 27 representation . Now, we are ready to calculate the potentials from 
Eq. (|4.ip with the profile function of Eq. (|4.3p . Fig. ^ plots the potentials for the fundamental, adjoint, and 27 
representations. The free parameters a, b, and / are chosen to be 0.05, 4, and 0.1, respectively. The potentials are 
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{14} = {3} @ {3} e {8} 

FIG. 6: Weight diagram of the adjoint (14-dimensional) representation of G(2) decomposed into its SU(3) subgroups. 

screened at large distances, as expected from the fact that G(2) has only a trivial center. In fact, the representations 
14 and 27 lead to a singlet after they combine with gluons that popped out of the vacuum. This happens when the 
distance between static sources is large enough and the potential between sources is big enough to create a pair of 
heavy gluons. 

{7} (S> {U} (E) {U} (E) {U} ^ {1} ® ... (5.8) 
{14}(8){14} = {1}® ... (5.9) 




{27} = {1} ® {3} ® {3} © {6} © {6} © {8} 

FIG. 7: Weight diagram of the representation 27 of G(2) decomposed into its SU(3) subgroups. 
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FIG. 8: Potentials between two G(2) static sources in the fundamental, 14 and 27 dimensional representations. The potentials 
are screened at large distances in agreement with the fact that gluons of the G(2) gauge group can couple to the initial sources 
and a singlet is produced. For each representation, a linear potential is observed at intermediate distances. 



{27} (E) {14} = {7} ® {14} ® {27} < 



(5.10) 



One of the interesting points of the figure is that the sources in the fundamental representation are screened at higher 
energy than the adjoint sources. This can be explained by the above equations which indicate that a fundamental 
quark interacting with three gluons leads to a singlet, while adjoint sources are screened by interacting with one 
gluon. Therefore, screening happens at higher energies for the fundamental representation compared to the adjoint 
representation. The sources of the 27-dimensional representation do not give a singlet after combining with a pair of 
gluons. Screening happens if the adjoint sources produced from Eq. (|5.10p interact with the adjoint sources again 
and a singlet is produced according to Eq. (|5.9|) . 

It is observed from Fig.® that the potentials are concave at intermediate distances. It is worse for representation 
14. To get a nonconcave potential and to be able to calculate the string tensions more accurately, we use the method 
of Ref. llj which considers some random fluctuations for the vortex profiles. Figure ([9]) shows the potentials at 



intermediate distances. A linear regime is observed for all three representations. However, the length of this area is 
very small. The Casimir scaling regime is expected to extend roughly from the onset of confinement to the onset of 
screening. String tensions are calculated for R between 5 and 12. From Fig. the string tension ratios are 



Kf 



= 1.48, 



while the Casimir ratios are 



Ci4 

c7 



= 2, 



K27 
Kf 



C27 



= 1.65, 



7 
3' 



(5.11) 



(5.12) 



The string tension ratios are qualitatively in rough agreement with Casimir ratios. The agreement between Casimir 
ratios and string tension ratios obtained from the lattice calculations is very good (^. However, we recall that the 
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FIG. 9: Potentials between two G(2) static sources in the fundamental, 14 and 27 dimensional representations at intermediate 
distances. The ratio of the string tension of each representation to the string tension of the fundamental representation is 
qualitatively in rough agreement with the Casimir ratio. 



thick center vortex model predictions for proportionality with Casimir scaling for SU(2), SU(3), and SU(4) have also 
been very rough [1, Q. The Casimir ratios of representations 27 to 14 is 1.16, which is in good agreement with the 
string tension ratios obtained from the model, about 1.12, in Eq. (|5.1ip . It is also possible to find the potentials 
between static sources in other higher representations of the G(2) gauge group, but finding the Cartan subalgebra 
from the weight diagrams would be a little bit harder. However, studying these three representations of the G(2) 
gauge group shows that the thick center vortex model, with the idea of the domain structure, works rather w ell even 
in the G(2) gauge group which has only a trivial center. 



VI. CONCLUSION 



Studying confinement in gauge groups without nontrivial center elements is an attractive subject. G(2) is one of 
these exceptional groups which has only a trivial center. Lattice calculations show some evidence for the confinement 
of quarks at intermediate distances. On the other hand, the thick center vortex model is a phenomenological model 
which gives potentials between quarks in gauge groups with nontrivial center elements. In this paper, to reproduce 
the lattice results, the idea of the vacuum domain structure has been used in the thick center vortex model. A flux 
has been assigned to the trivial center of the G(2), gauge group and the potential between two static G(2) quarks 
is obtained. The potential is screened at large distances, as expected from the thick center vortex model. In other 
words, the trivial center does not affect the Wilson loop, where a complete link is established between them at large 
distances. In this regime, the vortex (domain) is located completely inside the minimal area of the loop. A linear 
potential at intermediate distances is observed. The thickness of the domains and possibly the SU(3) subgroups of 
the G(2) gauge group are responsible for this linear behavior at intermediate distances. 

In addition, potentials between static sources of the adjoint and the 27 -dimensional representations are calculated by 
the model. In both cases, the potentials are linear at intermediate distances, and the string tensions are qualitatively 
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in rough agreement with Casimir scahng. 



VII. ACKNOWLEDGMENTS 



We would hke to thank Manfried Faber and Stefan Olejnik for very helpful discussions. We are grateful to the 
research council of the University of Tehran for supporting this study. 



[1] P. D. Collins, E. J. Squires, Regge poles in particle physics, 1968, Springer, New York. 

[2] G. 't Hooft, Nucl. Phys. B 138, 1978, p. 1; G. Mack, in Recent development in gauge Theories, edited by G. 't Hooft et 
al. (Plenum, New York, 1980); H. B. Nielson, P. Olcscn, Nucl. Phys. B 160, p. 380, 1979; J. Ambj0rn, P. Olcscn, Nucl. 
Phys. B 170, p. 60 and p. 265, 1980; J. Cornwall, Nucl. Phys. B 157, p. 392, 1979; R. Feynman, Nucl. Phys. B 188, p. 
479, 1981. 

[3] M. Faber, J. Greensite, S. Olejnik, Phys. Rev. D 57, p. 2603, 1998. 

[4] G. S. Bali, Casimir scahng of SU(3) static potentials, Phys. Rev. D 62, p. 114503, 2000. 

[5] S. Doldar, Phys. Rev. D 62, p. 034509, 2000. 

[6] L. Liptak, S. Olejnik, Phys. Rev. D 78, p. 074501, 2008; Bjorn H. Wellegchausen, A. Wipf, C. Wozar, Phys. Rev. D 83, p. 

016001, 2011; Bjorn H. Wellegchausen, A. Wipf, C. Wozar, Phys. Rev. D 83, p. 114502, 2011. 
[7] C. Bernard, Phys. Lett. B 108, p. 431; Nucl. Phys. B 219, p. 341, 1983; J. Ambj0rn, P. Olesen, C. Peterson, Nucl. Phys. 

B 240, p. 189, 1984; C. Michael, Nucl. Phys. B 259, p. 58, 1985; G. Pouhs, H. Trottier, Phys. Lett. B 400, p. 358, 1997. 
[8] S. Deldar, JHEP OlQl, p. 013, 2001, S. Deldar, S. Rafibakhsh, Phys. Rev. D 76, p. 094508, 2007. 
[9] J. Greensite, K. Langfeld, S. Olejnik, H. Reinhardt, T. Tok, Phys. Rev. D 75, p. 034501, 2007. 
[10] K. Holland, P. Minkowski, M. Pope, U. J. Wiese, Nucl. Phys. B 668, p. 207, 2003, M. Pcpc, U. J. Wiese, Nucl. Phys. B 

768, p. 21, 2007; Bjorn H. Wellegchausen, A. Wipf, C. Wozar, Phys. Rev. D 83, p. 016001, 2011. 
[11] S. Deldar, S. Rafibakhsh, Phys. Rev. D 81, p. 054501, 2010. 



